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Abstract We describe using Ap and F 10.7 as a geomagnetic-precursor pair to predict the 
amplitude of Solar Cycle 24. The precursor is created by using F 10.7 to remove the direct 
solar- activity component of Ap. Four peaks are seen in the precursor function during the 
decline of Solar Cycle 23. A recurrence index that is generated by a local correlation of 
Ap is then used to determine which peak is the correct precursor. The earliest peak is the 
most prominent but coincides with high levels of non-recurrent solar activity associated 
with the intense solar activity of October and November 2003. The second and third peaks 
coincide with some recurrent activity on the Sun and show that a weak cycle precursor 
closely following a period of strong solar activity may be difficult to resolve. A fourth peak, 
which appears in early 2008 and has recurrent activity similar to precursors of earlier solar 
cycles, appears to be the “true” precursor peak for Solar Cycle 24 and predicts the smallest 
amplitude for Solar Cycle 24. To determine the timing of peak activity it is noted that the 
average time between the precursor peak and the following maximum is ^ 6.4 years. Hence, 
Solar Cycle 24 would peak during 2014. Several effects contribute to the smaller prediction 
when compared with other geomagnetic-precursor predictions. During Solar Cycle 23 the 
correlation between sunspot number and F 10.7 shows that F 10.7 is higher than the equivalent 
sunspot number over most of the cycle, implying that the sunspot number underestimates the 
solar- activity component described by F10.7. During 2003 the correlation between aa and Ap 
shows that aa is 10 % higher than the value predicted from Ap, leading to an overestimate of 
the aa precursor for that year. However, the most important difference is the lack of recurrent 
activity in the first three peaks and the presence of significant recurrent activity in the fourth. 
While the prediction is for an amplitude of Solar Cycle 24 of 65 d= 20 in smoothed sunspot 
number, a below-average amplitude for Solar Cycle 24, with maximum at 2014.5 d= 0.5, we 
conclude that Solar Cycle 24 will be no stronger than average and could be much weaker 
than average. 
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1. Introduction 

Predicting solar activity is an active area of research. Our knowledge of the Sun and solar 
activity will be validated by testing such predictions, from the seconds that anticipate a flare 
to the decades needed for planning satellite missions. Numerous authors have advocated 
using geomagnetic precursors to predict the level of activity in an upcoming solar cycle. 
The value of a heavily smoothed geomagnetic index, such as aa or Kp, at solar minimum 
correlates with the strength of the next cycle (Ohl, 1966). To create predictions before so- 
lar minimum, geomagnetic-precursor pairs are created by removing the direct solar-activity 
component from the geomagnetic index, usually by a correlation fit with a solar- activity 
index. By removing the solar-activity component, the predictive power of the precursor is 
moved to an earlier phase of the solar cycle. Sargent (1978) and Ohl and Ohl (1979) devel- 
oped geomagnetic-precursor pairs by removing the geomagnetic activity directly correlated 
with solar activity from the geomagnetic activity index to produce predictions of Solar Cy- 
cle 21. Feynman (1982) suggested that the geomagnetic-activity index aa could be separated 
into one component proportional to the solar- activity cycle and an “interplanetary” compo- 
nent that resembles the sunspot cycle, but is delayed by several years after solar maximum. 
The interplanetary component often peaks just before solar minimum and was shown to be 
a better indicator for the amplitude of the following cycle by Hathaway, Wilson, and Re- 
ichman (1999). However, Brown (1992) showed that removing the solar- activity component 
did not produce better predictions of Solar Cycle 22. 

The frequency of solar flares and coronal mass ejections (CMEs), the main drivers of 
the solar-cycle component of geomagnetic activity, rises and falls with R z or F 10.7 (Feyn- 
man, 1982); updated through Solar Cycle 23 by Pesnell (2008a). Geomagnetic activity is 
also produced by the high-speed solar- wind streams flowing out of the coronal holes that 
develop out of phase with the sunspot cycle (Luhmann et al., 2002). During the decline to 
sunspot minimum, the low-latitude extensions of these coronal holes produce long-lasting 
high-speed streams that drive recurrent geomagnetic activity. The magnetic fields that give 
rise to these coronal holes may be the actual precursor of the upcoming sunspot cycle, al- 
though the causal relationship is not completely understood (Wang and Sheeley, 2009). 

Geomagnetic precursors were an important part of the consensus prediction for Solar 
Cycle 23 (Joselyn et al., 1997). However, all of the predictions for Solar Cycle 23 from 
geomagnetic precursor methods as evaluated by Hathaway, Wilson, and Reichman (1999) 
indicated a larger amplitude than was observed; although the actual amplitude was within 
or just outside their 2 o error estimates. The Solar Cycle 24 Prediction Panel also relied on 
geomagnetic precursors as an important component of their consensus prediction (Biesecker 
and the Solar Cycle 24 Prediction Panel, 2007). 

Predictions of the amplitude of Solar Cycle 24 using geomagnetic precursors cover a 
wide range of values (Pesnell, 2008b, 2012). One geomagnetic-precursor pair prediction 
combines aa and R z to predict a large amplitude for Solar Cycle 24 of 160 ± 25 for the 
maximum of the smoothed sunspot number (Hathaway and Wilson, 2006), although that 
value was later revised to below-average. A geomagnetic precursor using a time- shifted 
Kp predicts an amplitude of 96 (Kryachko and Nusinov, 2008). The trend of solar activity 
predicted by the geomagnetic-precursor pair of Hathaway and Wilson (2006) disagrees by a 
large amount with the predictions based on the solar polar field (Schatten, 2005; Svalgaard, 
Cliver, and Kamide, 2005). Reducing or understanding these discrepancies is necessary to 
give confidence to our ability to predict the amplitude of an upcoming solar cycle. 

Several issues could be present in the data used to construct a geomagnetic precur- 
sor of solar activity that would lead to an incorrect prediction. One is the adequacy of a 
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geomagnetic-activity index in representing that part of geomagnetic activity that could be 
a precursor. Another is how accurately the sunspot number [R z \ represents solar activity. 
Solar Cycles 22 and 23 both had secondary peaks in F 10.7 and a number of active regions 
in the southern hemisphere that were not as well represented in R z . We will show that the 
residual in a correlation fit between R z and F 10.7 shows that R z has been trending below 
F 10.7 since about 1996. Finally, we consider the effect of varying the width of the averag- 
ing filter and whether a weak-cycle precursor would be resolved when the filter width also 
includes a period of extremely high solar activity. 

In this article we examine how Ap and F 10.7 can be combined into a geomagnetic- 
precursor pair to forecast the level of activity of Solar Cycle 24. The predicted amplitude 
is of an averaged sunspot cycle, such as described by Hathaway, Wilson, and Reichman 
(1994), which can be described by a single number. Details such as Gnevishev’s gap, where 
two distinct peaks are observed, cannot be anticipated with this technique. The sensitivity 
of the predicted amplitude to the width of the averaging filter and association with recur- 
rent activity is explored. Next we show that the precursor should be evaluated at the time of 
peak recurrent activity near solar minimum. If this condition is not satisfied, the predicted 
amplitude will always be overestimated. The data series are described in the Appendix. 


2 . Ap and F10.7 as a Geomagnetic Precursor Pair 

Although the aa-R z precursor pair has been used by several authors to predict solar activity, 
we have six reasons to use Ap-Fio.7. 

2.1. Changes in R z Calibration 

One reason is the long-term reliability of R z . Svalgaard (2012) claims that R z was too 
small in the nineteenth century. If a long-term index has a bias toward being too large in 
the current epoch, then predictions using that index will show a similar bias. By using the 
well-calibrated F 10.7 as the solar- activity index we reduce these biases. 

2.2. Relative Behavior of F 10.7 and R z 

Another reason is the change in the relative behavior of F 10.7 and R z in Solar Cycle 23. 
If we assume an error of five in both F 10.7 and R z , there is a fairly linear correlation re- 
lationship between these two indices, given by R z = —63.8 + 1.07 F10.7. The residuals in 
this correlation fit [A R z \ contain two major contributions (shown by the blue line in Fig- 
ure 1). One is a trend for A R z to decrease from positive to negative during the period of 
time that F 10.7 is available. The other is an oscillatory component that was associated with 
other solar- activity observations until Solar Cycle 23. This residual shows that R z tends to 
be underestimated by F 10.7 in the ascending phase of the solar cycle and overestimated in 
the descending phase of the solar cycle. But A R z has remained negative (R z is smaller than 
predicted by the correlation fit with F10.7) since the solar minimum in 1996. If R z is too 
small in the current epoch, the solar- activity component of geomagnetic activity would be 
incompletely subtracted from aa, changing the magnitude of a precursor peak in that phase 
of the solar cycle. A higher value of the solar component in a precursor pair will tend to 
reduce the predicted level of solar activity. 

We have three reasons to interpret this trend as a decrease in R z rather than an increase in 
F 10 . 7 . Svalgaard (2009) investigated how F 10.7 tracks other radio-irradiance measurements, 
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Figure 1 The residual of the fit 
between Rz and F10.7 (straight 
annual averages.) The residual is 
shown as the blue solid line and 
the North-South asymmetry in 
the number of active regions is 
shown as a solid red line. The 
tendency of Rz to be smaller 
than F10.7 appears to be 
increasing with time. 
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primarily from the Japanese radio observatories. He concluded that the calibration of F 10.7 
can be traced through the entire time series and is not drifting to higher values in the current 
epoch. Tapping and Valdes (201 1) showed that several solar- activity indices, especially F 10.7 
and Rz, are well correlated at some times and less so at others. They associated these dif- 
ferences with the evolution of solar activity in the different layers of the solar atmosphere. 
Finally, Livingston, Penn, and Svalgaard (2012) also argued that the relationship between 
F 10.7 and R z appears to be drifting. 

A possible explanation for some of the variation is that F 10.7 tends to track the number 
of active regions in the Sun’s southern hemisphere during the second half of Solar Cycles 
22 and 23 while Rz does not. This is illustrated by the red curve in Figure 1, where the 
difference in the number of active regions in both hemispheres of the Sun (the N-S asym- 
metry), derived from the active regions reported for each year since 1950 by the National 
Geophysical Data Center (NGDC) website, is plotted. The number of active regions in both 
hemispheres of the Sun were derived from the active regions reported for each year since 
1965 on the NGDC website. Before 1965, the data values were obtained from the NGDC 
as yearly files containing the measured active-region locations and areas measured by the 
Greenwich Observatory. The overlap years were averaged. 

The residual A R z tracks the N-S asymmetry until 1996, after which it becomes in- 
creasingly negative. The residual had not returned to the =b 5 band at the end of 2008. This 
behavior is unlike the previous three cycles, where A R z tended to track the N-S asymmetry. 

A similar asymmetry in the location of flares, described as the relative difference of the 
northern and southern hemispheres, is described by Temmer et al. (2001). The asymmetry 
described here does have a definite variation with solar activity in the previous three cycles; 
early in the cycle there is an excess of northern active regions and later there is an excess of 
active regions in the southern hemisphere. A similar pattern is seen in the N-S asymmetry 
for the measured spot numbers since 1874 (Li, Gao, and Zhan, 2009). However, the south- 
ern hemisphere led the northern before Solar Cycle 19, while the North has led the South 
since Solar Cycle 20. Although the phasing is not perfect, fluctuations in the magnetic- 
field strength in sunspots shows a similar solar-cycle variation (Pevtsov et al ., 2013; 
McIntosh et al ., 2013). 

2.3. Relative Behavior of aa and Ap 

Next, aa is an antipodal index, where several effects are canceled by combining stations 
on diametrically opposite sides of the Earth. Over time the station pairs have changed and 
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offsets in aa have been suggested to maintain an accurate long-term index (Svalgaard, Cliver, 
and Le Sager, 2004). Ap is an average over many stations and is much less sensitive to 
changes in their locations. The residual in a correlation fit between aa and Ap shows a 
tendency for aa to be greater than the correlation fit during the decline of Solar Cycle 23, but 
the magnitude of that bias is very sensitive to whether or not the suggested offsets for aa are 
used. A lower value of the geomagnetic component in a precursor pair will tend to reduce the 
predicted level of solar activity. This effect is weaker than the A R z trend discussed above. 

2.4. Time Coverage 

A third effect is the importance of the temporal coverage, or number of previous cycles, 
of the precursor pair. Hathaway and Wilson (2006) show a precursor from about 1870 us- 
ing aa and Rz . However, their precursor has a more structured temporal dependence before 
1930 than after (there are 12 maxima in the precursor during the seven solar cycles between 
1870-1930; nine maxima in the eight solar cycles between 1930-2000). Another observa- 
tion is that all of the cycles before Cycle 17 have multiple peaks in the precursor and only 
one has a multiple peak after Cycle 17 (^ 1960, between Cycles 19 and 20). This could be 
a statistical fluctuation but it could also represent a true systematic effect in the data. A dif- 
ferent precursor pair can be examined to see if these effects are important. Because both Ap 
and F 10.7 are available since 1950, this pair can be examined over most of the modern period 
of the aa-i?z precursor function. Using F 10.7 as the activity indicator limits the number of 
previous cycles to five but still gives a correlation that spans a range sufficient to predict the 
amplitude of Solar Cycle 24. 

2.5. Effects of Smoothing Filter 

A fourth reason is that the wide smoothing filter advocated by Hathaway and Wilson (2006) 
encompasses both peaks of F 10.7 during solar maximum in Solar Cycles 22 and 23. The pres- 
ence and size of precursor peaks should be examined under other filter widths to see what 
effect this has on the predicted amplitude of Cycle 24. A precursor of a weak cycle will be 
quite small and could be disguised by the smoothing process. This was motivated by the Ap 
figure on the NOAA/SEC Solar Cycle Progression website (www.sec.noaa.gov/SolarCycle/), 
which shows peaks in 2005, 2006, and 2008 that may be more appropriate geomagnetic pre- 
cursors for Solar Cycle 24. 

2.6. Effect of Magnetosphere-Solar Wind Coupling 

Finally, what is the effect of the seasonal changes in the coupling between the terrestrial 
magnetosphere and the solar- wind fluctuations that drive geomagnetic activity? The extraor- 
dinary solar activity in late 2003 coincided with a maximum in the semi-annual Russell- 
McPherron effect (Russell and McPherron, 1973; Crooker, Cliver, and Tsurutani, 1992), 
when the tilt of the Earth’s rotation axis allows the most effective coupling between the 
magnetosphere and the solar wind. Without a clear understanding of this effect, it should at 
least be included in the uncertainty of the predicted amplitude. 


3. Construction of Ap-Fio .7 Geomagnetic Precursor Pair 

We adapted the methods originally described by Feynman (1982) and modified by Hath- 
away and Wilson (2006), but with Ap and F 10.7 as the data streams. The construction of a 
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Figure 2 The averaged Ap and scaled aa indices from 1930 through 2008. The aa index was scaled by 
the correlation fit Ap = — 4.00 + 0.807 aa and is listed as aa r in the legend. A wide (FWHM = 730 days) 
Gaussian filter was used to average the data. Note that aa r and Ap track each other until about 1980. This 
is particularly important for the geomagnetic-precursor analysis. Values of scaled aa are consistently higher 
than Ap since 1980, with the peak levels remaining constant while the peak values of Ap show a decrease 
over the same interval. The run of Fjo .7 (divided by ten) from 1950 through 2008 is shown for orientation. 


Figure 3 Averaged Ap and 
F 10.7 indices. The straight line is 
described by 

ApsA = 3.4 + 6.0 x 10 -2 Fjo. 7 - 
The fit represents the 
solar-activity contribution of Ap. 
After subtracting this fit from Ap, 
we obtain the geomagnetic 
precursor. 



geomagnetic-precursor pair using Ap and F 10.7 as the predictor/predicted pair is illustrated 
in Figures 2 through 5. 

The data were smoothed by the Gaussian filter recommended by Hathaway and Wilson 
(2006) with a width parameter [<r WG ] that was varied from 81 to 310 days. The widest 
filter had a FWHM of 730 days, encompassing almost four years of data when the filter 
was truncated to span ±2a WG . The filter was truncated at ±620 days and divided by the 
normalization of the truncated filter. The truncation at ± 2 <t wg means that 4.5 % of the 
filter was removed. Changing the width to ± 3a WG did not change any of the conclusions 
discussed here, but removing data at the recent end of the time series restricts how close to 
the present the precursor can be evaluated. 
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Figure 4 Temporal dependence of the interplanetary contribution to geomagnetic activity (in blue) that 
was found by subtracting the solar-activity component from Ap, which is then identified as the precursor. 
Excursions below zero were ignored. Maxima in both Fp ).7 and the precursor after 1950 are indicated. Values 
of these maxima were interpolated and used to construct the prediction points in Figure 5. The unaltered run 
of Ap is shown as a black line and the solar-activity component of Ap (from Equation (1)) is shown as a red 
line labeled as Fiq .7 in the legend. 


Figure 5 Predictor correlation 
fit between Ap and F 10 . 7 . The 
points marked by “+”s were 
interpolated from Figure 4. Cycle 
numbers are written near the 
plotted symbols. The straight line 
shows the linear-prediction 
relationship, whose fit 
coefficients are listed in 
Equation (2). For Ap p = 0.1 we 
obtain a predicted F 10.7 
amplitude of 122. Converting to 
annual-averaged sunspot number, 
this gives Rz,aa — 65, well 
below an average cycle of 
115 ±40. 



The wide-Gaussian averaged data for <j wg = 310 days are shown in Figure 2. A fit to the 
lower envelope of a plot of Ap vs. F 10 . 7 , 

A Psa = 3.4 + 6.0 X 10 " 2 F 10 . 7 , (1) 

is shown in Figure 3. We identified this fit as the solar- activity contribution to Ap and sub- 
tracted its value from each Ap value, leaving the interplanetary contribution, which is shown 
as the solid blue line in Figure 4. Maxima of F 10.7 and the precursor, indicated by downward 
arrows in Figure 4, are read from the graph and plotted in a correlation diagram in Figure 5. 
Solar cycle numbers are written near the plotted symbols corresponding to that cycle pair. 

The temporal dependence of the precursor in Figure 4 is similar to that of Figure 2 in 
Hathaway and Wilson (2006). One difference is that this peak value of the precursor de- 
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Table 1 Predictions of Solar Cycle 24: all have estimated errors of ± 20. 

Precursor [Ap p ] 

F 10.7,AA 

r 24,AA 

Method 

6.7 

183 

130 

awe = 310-day precursor evaluated at peak 
in 2003 

4.8 

165 

115 

awe = 150-day precursor evaluated at peak 
in 2005 

4.0 

158 

105 

cry/Q = 310-day precursor evaluated at time 
of maximum in 2005 from the 
o\y G = 150 day curve 

1.0 

130 

75 

awe = 150-day precursor evaluated at peak 
in 2006 

0.1 

122 

65 

o\y G = 150-day precursor evaluated at peak 
in 2008 from the or^ G = 81 -day curve 

Average 

125 

70 

Includes the last two predictions 


creases in a roughly linear fashion during the last three solar minima, while the aa -R z pair 
does not. The height of the second peak in the double peak around 1963 in the aa-/?z pair 
has been reduced in the Ap-Fio .7 pair. This could represent the asymmetry in R z where F 10.7 
does not show the same slope in the second half of the solar cycle as R z (see the Appendix 
and Figure 1). 

The straight line in Figure 5 shows the correlation relationship when the precursor for 
Cycle 20 is 4.2 (the second, smaller peak). A linear fit to the points in Figure 5, allowing for 
errors in both coordinates, is 

Fio.7.»+i = 120(±9) + 8.5(±1) A PpX 2 = 4.5. ( 2 ) 

Figure 4 shows that the precursor value in 2003 for Solar Cycle 24 with <x WG = 310 days 
is Ap p = 6.7. Using the linear fit between the annual- averaged and wide-Gaussian F 10 . 7 , 
Fio.7,aa = (Fio.7,wg — 13.3) /0.898, followed by the conversion between annual- averaged 
F 10.7 and R z [R z ,aa = 1.07 Fio. 7 ,aa — 63.8], this corresponds to R z ,aa = 130, above the 
amplitude of an average cycle [^z,ave = 1 15 =b 40 (Pesnell, 2008b)] but below the predicted 
amplitude of Hathaway and Wilson (2006). This is the first prediction listed in Table 1. 
The change from R z to F 10.7 and from aa to Ap has reduced the predicted amplitude by 30 
sunspot numbers. By encompassing over 1200 days in the average, it is impossible to look 
for precursors at later times. We now examine the behavior of the precursor with the width 
of the Gaussian smoothing filter. 


4. Variation with the Width of the Gaussian Filter 

One complication of the geomagnetic-precursor pair forecast method is determining the 
location in time of the precursor peak. One reason why the precursor peak occurred in late 
2003 and not closer to solar minimum could be the width of the Gaussian filter used in the 
analysis and the extraordinary level of solar activity in late 2003. This width also determines 
the temporal cutoff of the data. A filter with width a WG cannot be used to study effects nearer 
than about 2<t wg from the end of the time series. 

Versions of the Ap-Fio.7 precursor function using different widths of the Gaussian filter 
are shown in Figure 6. The widths were set to a WG = 81, 150, and 310 days, with a cutoff at 
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Figure 6 Versions of the Ap-F107 precursor function using different widths of the Gaussian filter. The 
widths were set to <r WG = 81, 150, and 310 days, with a cutoff at 2 <xwg before the end of the time series 
(the days column.) As crwG decreases, the value at the peak in 2003 increases, while a second peak in 2005, 
invisible in the o^g = 310 curve appears and becomes better resolved as decreases to 150. The shortest 
filter shows four distinct peaks from 2003 to 2008. The horizontal dotted lines show the prediction for the 
ct\VG = 310- and 81 -day precursors, the former at its peak and the latter at the peak coinciding with recurrent 
activity. 


Figure 7 The variation of the 
height and date of the precursor 
peak in mid-2003 to early 2004. 
The Ap-Fio .7 precursor function 
was evaluated with different 
widths of the Gaussian filter. As 
cr WG increases, the value at the 
peak in 2003 decreases, while the 
date of the peak moves toward 
later times. 



(days) 


=L2a WG (the days column). For each filter width the predictor function in Equation (2) was 
evaluated for all values of the precursor. This means that the predicted amplitude of Solar 
Cycle 24 can be read off the plot, leaving the choosing of the correct point in time as the 
remaining issue. 

Four peaks are present between 2003 and 2008 in the precursor using the narrowest filter 
length [a WG = 81 day]. As shown in Figure 7, the location of the peak in 2003 and its value 
changes in a systematic way with the filter width used in the analysis. The value of the peak 
in 2003 increases by 40 % as the filter width decreases by almost 30 % (Figure 7). The 
timing of the peak also shifts, moving back in time as the filter width is decreased. This 
correlated behavior of the peak and timing is due to the filter broadening to include two 
peaks of activity, one in 2003 and the other in 2005. The peak in 2005 can be seen in the 
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a WG =150 days, but decreases to a shoulder in the a WG = 210 curve, and becomes invisible 
when <j wg is increased to 310. 

Comparing the sunspot number predicted using the a WG = 150 and 310 day precursor 
values in 2005 (Ap p = 4.8 and 4.0, respectively) in the correlation functions of Equation (2) 
gives estimated amplitudes of Rz , aa =115 and 105, the amplitude of an average cycle. This 
variation implies that an uncertainty of at least ten sunspot numbers should be assigned to 
the prediction on this basis alone. 

Whether the Russell-McPherron effect is causing a modulation of the precursor cannot 
be addressed by this analysis. The width of the <r WG = 310 day filter (FWHM = 730 days) 
is wider than the seasonal dependence. When the filter length is reduced to 81 days, the 
seasonal dependence becomes more visible, but a lack of strict periodicity in this effect 
means that it is not averaged out by the filtering and must still be understood. The estimated 
uncertainty caused by this effect is eight in the predicted values of F10.7 and R z . 

Combining the formal error in Equation (2) (d= 10), the uncertainty in the value of the 
precursor (d= 10), and the variation from the Russell-McPherron effect [d= 8] in a root-mean- 
square gives an error estimate of d= 16, which we rounded to a prediction uncertainty of d= 20. 
The uncertainty in the timing of the precursor peak is about 0.5 years based on the width of 
the peaks in the narrowest filters (Figure 6). 

This analysis has not resolved which is the appropriate peak to use as a precursor of 
Cycle 24. The timing of the precursor could be related to the level of recurrent geomagnetic 
activity, which we now examine using a recurrence index. 


5. Comparison with a Recurrence Index 

The accuracy of the geomagnetic precursor is ascribed to its being sensitive to recurrent 
activity on the Sun. Geomagnetic activity near solar maximum is primarily driven by solar 
activity related to sunspots and active regions. These drivers are randomly distributed in 
longitude and relatively short lived, so the geomagnetic activity between successive solar 
rotations is uncorrelated. Even long-lived active regions will produce flares and CMEs at 
random times as they traverse the visible disk of the Sun. By comparison, the long-lived 
coronal holes present on the Sun during the decline from solar maximum are relatively fixed 
in heliographic longitude and emit high-speed streams that can strike the Earth every 27 
days, driving periodic geomagnetic activity that can be used as a proxy for recurrent activity 
on the Sun. This can be represented with a recurrent index [RI] using an Ap that is similar 
to the one defined by Sargent (1985) using aa as the geomagnetic indicator. 

We define RI as the correlation in Ap from one Bartels rotation to the next: 


where (Ap)/ is the average of Ap over Bartels rotation i, 07 is the standard deviation from 
the average in that interval, and the time assigned is the beginning of rotation /, the middle 
of the combined intervals. (A Bartels rotation number is the number of 27-day solar-rotation 
periods since the first began on 8 February 1932.) RI and the Ap-Fio.7 precursor from 1950 
through 2008 are plotted in Figure 8. The level of recurrent activity is highest near solar 
minimum and lowest near solar maximum. 

Each of the precursor peaks prior to Solar Cycles 20 - 23 coincides with a local max- 
imum in RI exceeding 0.3. If recurrent activity were the source of the precursor peak in 
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Figure 8 The recurrence index [RI] defined by the correlation of Ap between one Bartel rotation and the 
next (Equation (3)). The curve was smoothed over 15 rotations. Shown as a blue dashed line is the Ap-Fio .7 
precursor function from Figure 4. Vertical blue dashed lines show the maxima of the precursor used to con- 
struct the correlation function in Figure 5. The vertical red dotted lines show the minima in 2006 and 2008. 
High values of RI (approaching and exceeding 0.5) are seen during solar minimum, while RI is close to (and 
sometimes can dip below) zero at solar maximum. Solar activity was dominated by several large active re- 
gions in the second half of 2003, which reduced RI to near zero. As we neared solar minimum in 2008, RI 
reached values higher than previously seen, but has since declined. A recurrence index derived from aa shows 
similar behavior. 


2003, a peak in RI should coincide with that precursor peak. The temporal variation of RI 
since 2000 and the Ap-Fio.7 geomagnetic-precursor pairs for different smoothing widths are 
plotted in Figure 9. A local minimum in the recurrent index is indeed present in late 2003, 
with RI approaching zero. If the filter width is increased to <r WG = 365 days, the peak moves 
into early 2004, where RI is increasing but remains lower than 0.2. Even the second possible 
precursor peak in 2005 (calculated with the narrower filters) is not coincident with a local 
maximum of recurrent activity. It is not until 2008 that unambiguous recurrent geomagnetic 
activity is present, and it is at this time that a fourth peak appears in the narrowest filters. 
This precursor is Ap p = 0.1 and gives Fio. 7,24 = 122 (Table 1). 


6. Conclusions 

Using Ap and F 10.7 as a geomagnetic-precursor pair gives a predicted amplitude for Solar 
Cycle 24 that is smaller than a similar analysis using aa and R z . Once the interplanetary 
component of Ap is derived and used as a precursor, we obtained an F 10.7 amplitude of 
122 =b 20 and a sunspot number of 65 d= 20 for Solar Cycle 24, with the peak of solar activity 
in May 2014. The combination of R z underestimating F 10.7 and aa overestimating Ap in 
the latter parts of Solar Cycle 23 leads to a prediction of a lower-amplitude cycle when Ap 
and F 10.7 are used as the precursor pair. Requiring significant recurrent activity to be present 
also delays the timing of the peak and allows the precursor to drop well below the 2003 
value. This prediction for a below-average Solar Cycle 24 agrees with the SOlar Dynamic 
Amplitude (SODA) Index prediction of Schatten (2005) and others that use the solar polar 
magnetic field as a precursor of future solar activity (Svalgaard, Cliver, and Kamide, 2005). 

We have examined how the precursor’s value and location varies with the width of the 
smoothing filter and how the presence of multiple peaks in the precursor can lead to misiden- 
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Figure 9 Versions of the Ap-Fio .7 geomagnetic precursor using widths of <r W G = 81, 150, and 310 days. 
Plotted as a dashed line is the recurrence index defined by Equation (3) smoothed over 15 rotations. The 
horizontal dotted line shows the zero level of RI. If recurrent activity were the source of the precursor peak in 
2003, a peak in the recurrent index should coincide with the precursor peak in 2003. A local minimum in the 
recurrent index is indeed present at that time, with the index approaching zero. Smoothing of the recurrent 
index over more rotations does not remove the dip in 2003. Even the second possible precursor peak in 2005 is 
not coincident with a local maximum of recurrent activity. The widest peaks represent averages over episodes 
of non-recurrent and recurrent activity. That tends to increase the precursor by emphasizing the extremely 
high levels of activity in late 2003 and smearing the later, narrower peaks. 

tifying the timing of the precursor peak. This ambiguity of multiple precursor peaks can be 
resolved by requiring that substantial recurrent activity coincides with the peak. This allows 
narrower filters to be applied during the smoothing step, reducing edge effects and revealing 
weaker features. 

Brown (1992) showed that predictions of the amplitude of Solar Cycle 22 using bivariate 
functions (such as the aa -R z and Ap-Fio.7 precursor pairs discussed here) were less accurate 
than those using only the geomagnetic index at solar minimum as the predictor function. 
This may be repeated in the amplitude predictions of Solar Cycle 24. However, by delaying 
the evaluation of the precursor function until almost solar minimum, the Ap-Fio.7 precursor 
is reduced to the minimum- Ap precursor. 

From a comparison with an index of recurrent geomagnetic activity, we showed that the 
peak of the precursor function in 2003 is associated with non-recurrent geomagnetic activity. 
In previous cycles the precursor peak coincided with a peak in the recurrence index defined 
in Equation (3). Two more precursor peaks, accompanied by some recurrent activity, were 
found in 2005 and 2006 using narrower filter widths. The behavior of the Solar Cycle 24 
precursor would then be similar to the one preceding Solar Cycle 20, another small cycle. 
That is not true of the precursor peak in 2003, which is dominated by activity with a very low 
value of RI. Our narrowest filter (<x WG = 81 days) has a peak in early 2008 with a value of 
approximately 0.1 that coincides with RI ~ 0.6. When used in Equation (2), this corresponds 
tO R24,AA = 65. 

Another precursor prediction uses the strength of the polar magnetic fields to estimate the 
amplitude of an upcoming solar cycle. Using this technique, Svalgaard, Cliver, and Kamide 
(2005) and Schatten (2005) predicted small amplitudes (75 =L 8 and 80 =L 20, respectively) 
based on the weak polar fields observed on the Sun during the decline of Cycle 23. This 
precursor is based on solar-dynamo theory, where the polar fields at solar minimum form 
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the seed for the next cycle’s toroidal fields and activity level. Predictions for Cycle 23 using 
the polar-field precursor also tended to be lower and in better agreement with the actual 
amplitude of Cycle 23 than the geomagnetic precursor. 

The earliest SODA Index prediction for Solar Cycle 24 was reported by Schatten (2005). 
SODA is designed to be a continuously updated prediction (Schatten and Pesnell, 1993). 
Looking at the SODA Index over the years between the first prediction and solar minimum 
gives an average prediction of 138 ± 2.5 for F 10.7 or 84 d= 2.5 for R z . This consistency in 
the prediction is in contrast to the change in the geomagnetic precursor prediction of 50 
sunspot numbers (40 %) during the same period of time. By relying on a particular time 
in the decline of the solar cycle, geomagnetic-precursor pairs tend to be a upper limit until 
much closer to solar minimum than is the polar field precursor. 

Solar polar magnetic fields are a direct measurement of one aspect of the solar dynamo, 
and it is reasonable to expect a precursor relationship between those fields and future solar 
activity. In addition, the solar polar field strength was slowly varying near a peak value in 
the recent solar minimum, lending credence to a prediction using that field strength. Geo- 
magnetic precursors are at least one layer of abstraction removed from the solar magnetic 
fields that presage the next cycle. Regardless of which solar magnetic fields are the source 
of the geomagnetic activity, the uncertainty of that coupling layer may be increasing. An 
understanding of this layer, which could drift in time, is needed to link the response of the 
Earth to the future behavior of the Sun. Geomagnetic activity also tends to decrease to very 
low levels at solar minimum, adding to the difficulty of using aa and Ap as a precursor. 

Whether F 10.7 or R z more accurately represents solar activity for use in geomagnetic 
precursors can only be assessed by knowing that they currently differ and monitoring their 
future behavior. This research has shown that having solar- activity indices that are accurate 
and cover long times is an important part of solar-cycle prediction. As predictions of the 
next solar cycle are made at ever-earlier times before solar minimum, the sensitivity of geo- 
magnetic precursors to filter length is quite troubling. A requirement of coincident recurrent 
activity may ensure that a proper prediction is made. 
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Appendix: Sources and Preparation of Data 

Predicting solar activity with a geomagnetic-precursor pair requires an index of geomagnetic 
activity and an index of solar activity. The solar-activity indices and two geomagnetic indices 
used in this article are described here. 

Rz 

The sunspot number [R z \ is often used as a long-term index of solar activity. It has been 
measured or derived for roughly 400 years and is the source of much of our knowledge of 
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the evolution of solar activity. Several versions of the sunspot-number record, such as the 
International, American, and group, are available; we used the International. Problems in 
the homogeneity of the R z time series are described by Cliver and Svalgaard (2007). 

Fio.7 

The solar spectral irradiance at a wavelength of 10.7 cm [F 10 . 7 ] is a well-known index of 
solar activity, which has been measured on a daily basis since 1947 (Covington, 1969; 
Tapping and Charrois, 1994). F 10.7 peaks at solar maximum and fades to around 60 x 
10 -22 W m -2 Hz -1 = 60 SFU at solar minimum. Because it has a latency of 24 hours, F 10.7 
is often used as a proxy for solar activity in models of the thermosphere (to understand satel- 
lite drag) and ionosphere (for radio- wave interference). Therefore, predictions of F 10.7 are 
very useful to operators of satellites and the space- weather community. 

aa 

The aa index is the average of magnetic activity measured at two antipodal observatories at 
an invariant magnetic latitude [A] of about =b 50° but located at diametrically opposite sites 
on the Earth (Mayaud, 1980). The antipodal average is used to approximately cancel the 
diurnal and annual variations in the aa index (Mayaud, 1972). 

Ap 

Another widely used geomagnetic index is the daily equivalent planetary amplitude or Ap 
index. It is constructed using amplitudes from 13 stations to approximately account for the 
local time and annual variations. The aa and Ap indices have a similar temporal dependence 
because they are both based on the three-hourly K indices measured at the stations in the 
same network of geomagnetic observatories (Menvielle and Berthelier, 1991), but Ap is 
derived from a larger network of stations. 

Active Region Data 

The number of active regions in both hemispheres of the Sun were derived from the ac- 
tive regions reported for each year since 1965 by the National Geophysical Data Center 
(NGDC) website. Before 1965 the data values were obtained from the NGDC as yearly files 
containing the measured active-region locations and areas from the Greenwich Observa- 
tory (ftp.ngdc.noaa. gov/STP/SOLAR_DATA/SUNSPOT_REGIONS/Greenwich/). The over- 
lap years were averaged. A summary file of annual active-region areas in the northern 
and southern hemispheres (ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SUNSPOT_REGIONS/ 
Greenwich/docs/ss_area.n_s) showed similar behavior. 
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